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Abstract. Ramanujan sums are exponential sums "with exponent defined over the 
irreducible fractions. Until no"w, they have been used to provide converging expansions 
to some arithmetical functions appearing in the context of number theory. In this 
paper, we provide an application of Ramanujan sum expansions to periodic, quasi- 
periodic and complex time series, as a vital alternative to the Fourier transform. The 
Ramanujan-Fourier spectrum of the Dow Jones index over 13 years and of the coronal 
index of solar activity over 69 years are taken as illustrative examples. Distinct long 
periods may be discriminated in place of the 1//" spectra of the Fourier transform. 



PACS numbers: 02.10.De, 05.45.Tp, 89.20-a 
1. Introduction 

Signal processing of complex time- varying series is becoming more and more fashionable 
in modern science and technology. Indices arising from the stock market, changing 
global climate, communication networks such as the Internet etc. are widely used tools 
for business managers or governmental representatives. There already exists a plethora 
of useful approaches for signal processing of complex data. The oldest and perhaps most 
widely used method is a Fourier analysis and its "fast" implementation: the fast Fourier 
transform (or FFT). Other complementary techniques such as wavelet transforms, 
fractal analysis and autoregressive moving average models (ARIMA) were developed 
with the aim of identifying useful patterns and statistics in otherwise seemingly random 
sequences p. 

Ramanujan sums are defined as power sums over primitive roots of unity. One can 
use an orthogonal property of these sums (closely related to the orthogonal property 
of trigonometric sums) to form convergent expansions of some arithmetical functions 
related to prime number theory [21 ISj. Following the ideas of Gadiyar and Padma 
[1], the first author proposed to expand the domain of application of Ramanujan sum 
analysis from number theory to arbitrary real time series and introduced the concept 
of a Ramanujan-Fourier transform [5j. This earlier work remained quite ambiguous 
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about the detection of isolated periods. Ramanujan sum expansions of divisor sums, 
sums of squares, and the Mangoldt function are well known. Surprisingly, the detection 
of a singly periodic signal by the Ramanujan sum analysis has not been considered 
before. But the Ramanujan- Fourier amplitude corresponding to a single cosine function 
of period q is extremely simple: as we shall see, the amplitude of the cosine function is 
simply scaled by the cosine of the delay and the inverse of the Euler totient function 
0(no). Similarly to the standard discrete Fourier transform, there are spurious signals 
of magnitude O(no/t), depending of the length t of the averaging spectrum. 

In the discrete Fourier transform, a sample to be analyzed is discretized into pieces 
of length 1/q and the expansion is performed over the g-th complex dimensional vectors 
of the orthogonal basis e'q\n) := exp (^^n), {p = 1, . . . , g). The orthogonal property 

reads eg^''(n)eg*'*(n) = q6{r, s), where 6{r, s) is the Kronecker symbol. The expansion 

of a time series is a{n) = ^pO'peq\n) with Fourier coefficients ctp = ^ Xln 
where the summation runs from to g — 1. In the Ramanujan- Fourier transform, the 
expansion a{n) = ^gCLqCqin) over the Ramanujan sums Cq{n) = Yl'p^q'\'^) (see Sec.|2]) 
involves the resolution - at every single scale from g = 1 to t — oo. The deep principle 
behind rests on a very intricate link between the properties of irreducible fractions | 
and prime numbers [2]-[S]- a result, one finds a much finer structure of time series, 
with a variety of novel features. 

The paper is organized as follows. In Sec. [2], we remind the reader with the 
arithmetical properties of Ramanujan sums, provide the definition of the Ramanujan- 
Fourier transform and examine the detection of a cosine signal. In Sec. El the use of the 
method is illustrated on the data from the stock market and solar activity. 

2. Ramanujan sums and the Ramanujan-Fourier transform 

Ramanujan sums are real sums defined as n-th powers of g-th primitive roots of the 
unity. 



where the summation runs through the p's that are coprime to g (hence the use of the 
symbol ""'), being first introduced in the context of number theory jS, [3] for obtaining 
convergent expansions of some arithmetical functions such as the relative sum of divisors 
a{n)/n of an integer number n, 



They are multiplicative when considered as a function of g for a fixed value of n, which 
can be used to prove an important relation 
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q=l 
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Cg{n) = /^'^^/^i)^^^' ^1 = 

In the above relation, the Euler totient function 0(g) is the number of positive integers 
less than q and coprime to it. The Mobius function, fi{n), vanishes if q contains a square 
in its (unique) prime number decomposition Hi "yi"' ili ^ prime number), and is equal 
to (—1)^ if q is the product of k distinct primes. One can readily checks the following 
orthogonal property 
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'^^Cr{n)cs{n) = 1 if r = s and ^^Cq(n) = q(f){q) otherwise. 

ra=l n=l 

For an arithmetical function a{n) possessing a Ramanuj an- Fourier expansion 



^aqCq{7 



a[n) 

with Ramanuj an- Fourier coefficients a^, one can write a Wiener-Khintchine formula, 
relating the autocorrelation function of a(n) and its Ramanuj an- Fourier power spectrum, 

^ oo 

limTv^oo^ 5Z "(^)«('^ + ^) = «?Cg(/i). 

n<N 9=1 

This relation was used for counting the number of prime pairs within a given interval 
[3]. A similar formula has been proposed for the convolution and cross-correlation [6]. 

Clearly, the Ramanuj an sum analysis of an arithmetical function looks like the 
Fourier signal processing of a time series a{n) at discrete time intervals n. This 
formal analogy was developed in [5] for the processing of time series with a rich low 
frequency spectrum [7]. Ramanujan signal processing was further developed in the 
context of quantum information theory [8] eventually leading to an original approach 
of quantum complementarity [9]. The Ramanuj an- Fourier transform was also used for 
processing time series of the shear component of the wind at airports [10], the structure 
of amino-acid sequences [11] and in relation to the fast Fourier transform [12]. All these 
applications make use of the property that for arithmetical functions possessing a mean 
value 



A^{x) = limt^ooy ^ 
one can write the inversion formula 



t 
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aq = -jj-:;Aj,{a{n)cq{n)). 
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Figure 1. Ramanujan- Fourier spectrum for three different cosine functions of periods 
no = 10, 14 and 30, computed from a sample of length t = 100. The amplitude at 
q = no equals l/^(no) = 1/4, 1/6 and 1/8, respectively. 
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Figure 2. Ramanujan- Fourier spectrum for the cosine function of period uq = 38, 
with delays 5 = 0, n/2 and tt, and sample length t = 100. The amphtudes of the peaks 
for 5 = and tt is 1/6. 




Figure 3. Ramanujan-Fourier spectrum for the cosine function of period no = 18, 
with sample lengths t = 100 and 500. One clearly observes the compression of the 
lines when t increases. 

Ramanujan-Fourier transform of a cosine function 

Let us consider now the Ramanujan signal processing of a periodic (cosine) function of 
period uq 

Tl 

a{n) = ao cos(27r h 6). 
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The Ramanuj an- Fourier coefficients read 

aoexp{i5) ^J-^ p 1 

Uq = hmt-^^——— > ^ > ^ exp 2z7m(- H ) + c.c. 

20(g)t q no 

in whicli tlie order of summations is reversed and "c.c." stands for tlie complex conjugate. 
Assume first that t is a muhiple of the period Uq. Then the n-th summation is zero 
unless /c = I + ^ is a positive integer. For instance, can be non-zero if q = Uq under 
the conditions that Uq divides p + 1 and k = 1, i.e. p = Uq — 1 (otherwise p > q, which 
is outside the range of summation of the p sum) . One then gets 

aoexp(i5)^ n 
«no = limt^oo^TTT — > ^exp(2z7r— ) + c.c. 

The n-th summation equals zero unless Uq divides n, otherwise it equals to As a 
result, the amplitude of the no-th line reflects the amplitude of the periodic signal as 

■cos(oj. 



"° ^(no) 

In general, t is not a multiple of uq so that there exists an extra contribution to the 
amphtude, of order of magnitude O(^). As long as the period no is much smaller than 
the length t of the sample, i. e. no ^ one observes a single line at no; otherwise bursts 
of non-zero amplitudes emerge in the vicinity of the lines kriQ — see Figs. 1-3. 

Thus, there arc two significant differences when compared to a period analysis by 
the standard discrete Fourier transform. First, the amplitude of the line at the period 
no is scaled by a factor of (f)(no). Second, the Ramanuj an- Fourier analysis is sensitive 
to the delay 5. The latter feature may, at first sight, seem as a drawback since some 
period of the signal to be analyzed may be hidden by the dephasing effect. One method 
to circumvent this difficulty is to average the spectra corresponding to several shifted 
samples of the signal. 



Ramanuj an- Fourier transform of a period modulated cosine function 

Let us now apply the approach to a period modulated cosine function. We intentionally 
select a period modulation with a large index (equal to 1). The selected modulation is 

71 

no = no[l + sin(27r — )], 

Til 

with no = 10 and ni = 14. The sample length is t = 2000. Due to a high modulation 
index, the FFT analysis (shown in Fig. 4) does not easily allow to recover the constituent 
integer periods 10 and 14. In contrast, the Ramanujan sum analysis is very powerful in 
this context. From Fig. 5 one clearly identifies (positive) large amphtudes at the periods 
10, 12, 2 X 12 and LCM(10, 12) = 70 (LCM being the least common multiple). Thus, 
for an input signal of the period no and period modulation ni, the FFT exhausts all 
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lines at Ihq + mni (/ and m integers), eventually leading to a continuous spectrum in 
the limit of incommensurate periods uq and rii. In contrast, the Ramanuj an- Fourier 
transform is straightforward in identifying the input modulation. 
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Figure 4. FFT spectrum of a period modulated cosine: no = 10 and ni = 14. 
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Figure 5. Ramanuj an- Fourier transform of a period modulated cosine function: 
no — 10 and n\ = 14. 



3. Ramanujan sums analysis of some complex systems 

As a nice illustration of the above-outlined properties the Ramanujan sum analysis, we 
shall analyze a couple of complex time series taken form the stock market and solar 
activity. 

The Dow Jones index of the stock market 

The first time series deals with the Dow Jones index and has been downloaded from 
http : // www.optiontradingtips.com/resources/historical - data/dow-jones30.html. 
Fig. 6 depicts the evolution of Dow Jones 30 Industrials stock price over about 13 
years. The power spectral density of the prices (Fig. 8) approximately follows a 1//^ 
law versus the Fourier frequency / = n~^, compatible with a Brownian- motion-based 
model [131 [H]- The Ramanuj an- Fourier analysis shown in Fig. 7 yields a more detailed 
structure with many (positive or negative) peaks centered at well-identified frequencies. 
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There exists a sensitivity of the amphtude of the peaks (not shown) on the number t 
of data, but the position of the peaks, as well as their statistics, is not dependent on 
t. Since both spectra in Figs. 6 and 7 are given in a logarithmic time-scale, it follows 
that the Ramanujan sum analysis provides a clear advantage over the standard Fourier 
analysis in offering a rich and structured signature. Here, we shall not delve any further 
into the origin of this structure, which will be a topic of a separate paper. 
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Figure 6. The Dow Jones 30 Industrials from 03.01.1995 to 30.05.2008. 
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Figure 7. The Ramanujan sums analysis of the Dow Jones 30 Industrials. Periods in 
the range n = 100 to n = 1000 were selected. One clearly sees (positive and negative) 
peaks centered about non-equally-spaced periods. 




Figure 8. The FFT analysis of the Dow Jones 30 Industrials in a log-log scale. 

The standard 1//^ dependence of the power spectrum, which is characteristic of a 
Brownian-like motion, is clearly visible. 
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The coronal index of the solar activity 

The second time series has been picked up from 

http : // www.ngdc.noaa.gov/stp/SOLAR/ftpsolarcorona.htmlT^index. 
It represents the Green Line (FeXIV 530.3 nm) Coronal Index of solar activity from 
1939 to 2008. One easily recognizes from Figs. 9 and 10 that the coronal index 
is approximately periodic, with a period about 10 years. The whole FFT spectrum 
shown in Fig 11 exhibits a 1// dependence characteristic of many physical, biological, 
arithmetical [TJ [5l [T5j and other complex systems [15] . 
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Figure 9. The temporal variation of the Coronal Index; a 10-year period is clearly 
visible. 
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Figure 10. The RFT of the Coronal Index. The 10-year period is clearly identified; 
other longer periods of a smaller amplitude are observed as well, besides the harmonics 
priQ (p integer). 
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Figure 11. The FFT of the Coronal Index. One observes a 1// dependence of the 
power spectrum; the 10-year period is quite hard to identify. 

Perspectives and Conclusion 

It is a widely shared belief that 1//" noises are so random that non-statistical models 
of them are currently out of reach. A counterexample to this belief can be found in [7j, 
in which an arithmetical approach to 1// noise was suggested. In the present paper, we 
offer another perspective by analyzing the data from an arithmetical magnifying glass 
built on Ramanujan sums. The Ramanuj an- Fourier transform is able to extract quasi- 
periodic features which are characteristic of number theoretical functions [2]-[5], as well 
as fine periodic features that the standard Fourier transform may hide. A Ramanujan 
sums analysis is a multi-scale prism with scales related to each other by the properties of 
irreducible fractions. It is particularly well-suited for analyzing rich time series showing 
a 1//" (0 < a < 2) FFT dependence. We selected two specific complex systems to 
illustrate the power of this new method: the data from the stock market (for which 
the price index FFT follows a 1/ /^-law) and those from solar cycle activity (for which 
the coronal index follows a l//-law). A more detailed examination of the latter will be 
given in a separate paper. 
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